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P-stars are compact stars made of up and down quarks in /3-equilibrium with 
J> ■ electrons in a chromomagnetic condensate. P-stars are able to account for com- 

pact stars as well as stars with radius comparable with canonical neutron stars. 
We compare p-stars with different available observations. Our results indicate 
^ ■ that p-stars are able to reproduce in a natural manner several observations from 

^ ■ isolated and binary pulsars. 
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■ 1. INTRODUCTION 



^ . In few years since their discovery (iHewish et al. 19681 ) pulsars have been identified with 



rotating neutron stars, first p redicted theoretically by Baade & Zwicky (1934a,b,c), endowed 



with a strong magnetic field (jPacini 19681 ; iGold 19681 ) . 



It is widely believed that there are no alternative models able to provide as satisfactory 
an explanation for the wide variety of pulsar phenomena as those built around the rotating 
neutron star model. Nevertheless, recently we have proposed (Cea 2004a,b) a new class of 
compact stars, p-stars, which is challenging the standard paradigm. 

P-stars are compact stars made of up and down quarks in /^-equilibrium with electrons 
in a chromomagnetic condensate. In our previous studies (Cea 2004a,b), we find that p-stars 
are able to account for compact stars, stars with radius comparable with canonical neutron 
stars, as well as super massive compact objects. Moreover, we show that p-stars once formed 
are absolutely stable. The cooling curves of p-stars compare rather well with observations. 
We also suggested that p-matter produced at the cosmological deconfinement phase transi- 
tion could be a viable candidate for baryonic cold dark matter. 
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In addition, in Cea (2006) we discuss p-stars endowed with super strong dipolar mag- 
netic field. We find that soft gamma-ray repeaters and anomalous X-ray pulsars can be 
understood within our theory. In particular, we point out that within our p-star model there 
is a quite natural mechanism to account for the generation of dipolar surface magnetic fields 
up to 10 16 Gauss. We succeed in obtaining a well defined criterion to distinguish rotation 
powered pulsars from magnetic powered pulsars. We show that glitches, that in our theory 
are triggered by magnetic dissipative effects in the inner core, explain both the quiescent 
emission and burst activity in soft gamma-ray repeaters and anomalous X-ray pulsars. We 
are able to account for braking and normal glitches observed in soft gamma-ray repeaters 
and anomalous X-ray pulsars. We discuss and explain the observed anti correlation between 
hardness ratio and burst intensity. Within our p-star theory we are able to account quanti- 
tatively for light curves from both gamma-ray repeaters and anomalous X-ray pulsars. 

It is worthwhile to briefly discuss the theoretical foundation of our proposal. Indeed, 
quite recently, the QCD vacuum was probed by means of an external const ant abelian chro- 
momagnetic field on the lattice (ICea &: Cosmai 20031 ; ICea fc Cosmai 20051 ). We found that 
by increasing the strength of the applied external field the deconfinement temperature de- 
creases towards zero. In other words, there is a critical field gH c such that for gH > gH c 
(g is the color gauge coupling and H is the strength of the chromomagnetic field directed 
along the third direction in color space) the gauge system is in the deconfined phase (the 
color Meissner effect). As a consequence, we see that there is an intimate connection between 
abelian chromomagnetic fields and color confinement. The existence of a critical chromomag- 
netic field is compatible with the QCD vacumm behaving like a disordered chromomagnetic 
condensate, for strong enough chromomagnetic field strengths enforce long range colour or- 
der thereby destroying confinement. In a seminal paper, R. P. Feynman (1981) argued that 
QCD in two spatial dimensions the confining vacuum at large distances is a chromomag- 
netic condensate disordered by the gauge invariance. The confinement of colors comes from 
the existence of a mass gap and the absence of colour long range order. As a matter of 
fact, we were able to show that, indeed, the QCD vacuum in three spatial dimensions does 
display at large distance a mass gap and no colour long range order and that the vacuum 
displays the color Meissner effect. We find that chiral symmetry breaking for quarks in 
the fundamental color representation is an inevitable consequence of the disordered chro- 
momagnetic condensate. Finally, we argue that the deconfined vacuum behaves more like 
a correlated liquid than an ideal gas, in accordance with recent results from heavy ion ex- 
periments. This, in turns, leads us to the conclusion that the deconfined QCD vacuum 
is characterised by long-range chromomagnetic correlations and whence p-matter, namely 
almost massless up and down quarks immersed in a chromomagnetic condensate, must be 
quite close to the true QCD deconfined state. In absence of strong gravity effects, the chro- 
momagnetic condensate in p-matter is of order of the critical field strength, which turns out 
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to be 4gB~c ^ 1.0 Gev flCea fc Cosmai 20031 : ICea fc Cosmai 2005T ). 

Note that we are using natural units where H — c — 1. In this units g is dimension- 
less and \f gH has dimension of energy. To switch to cgs units it suffices to replace gH with 
he gH. 

On the other hand, when including the effects of gravity it turns out t hat p-matte r gives 
rise to compact stars (p-stars), which are more stable than neutron stars (ICea 2004al ) what- 
ever the value of \fgH- Therefore, we have in general that in p-stars the chromomagnetic 
condensate must satisfy the constraint: 



gH < Jg~W c ~ 1.0 Gev . (1) 

We see, then, that the chromomagnetic condensate in p-stars acts like a dynamical effective 
bag constant which can be varied according to Eq. (TjQ). The actual value of the chromomag- 
netic condensate depends only on the central density of the compact star. In this way, p-stars 
are able to overcome the gravitational collapse irrespective to the stellar mass. Indeed, as 
discussed in Cea (2004a), our p-stars do not admit the existence of an upper limit to the 
mass of a completely degenerate configuration. In other words, our peculiar equation of state 
of degenerate up and down quarks in a chromomagneticn condensate allows the existence of 
finite equilibrium states for stars of arbitrary mass. To s ee this, we note that on dim ensional 
ground from the Tolman-Oppenheimer-Volkov equation (IShapiro fc Teukolsky 19831 ) we get: 
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1 



f(Sc 



R 



where e r 



G 3/2 gII •> v cj G 1 ' 2 gH 
e c / {gH) 2 and e c is the stellar central density. As a consequence we obtain: 
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(3) 



For reader's convenience, it is worthwhile to rewrite the previous equations in cgs units. 
We get: 
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where p c 



G s / 2 gH J ^ CJ G 1 / 2 gH 
(h c 3 p c )/(gH) 2 and p c is the stellar central density in gr/cm 3 , 
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(4) 



(5) 



From Eqs. (j^J) and we see that by decreasing the strength of the chromomagnetic 
condensate we increa se the mass and radius of the star, while the ratio 2 R M depends only 
on e c . It turns out (ICea 2004al ) that the function h(x) defined in Eq. ([3]) is less than 1 
for any allowed values of e c . By varying the chromomagnetic condensate strength in the 
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allowed range, the M-R p-star curves span the allowed region in the M-R plane. Indeed, a 
clear observational evidence in favour of p-stars would be the detection of pulsars with mass 
above 3 Mq. 

In the present paper we focus on few selected topics with the aim to furnish further com- 
pelling evidences in support of our proposal. In section [2] we compare recent determination 
of masses and radii of isolated and binary pulsars with our p-star model. In section [3j after 
a brief review of cooling in p-stars, we compare our theoretical cooling curves with several 
observational data. Finally, we draw our conclusions in section HI 



2. MASSES and RADII 

In a recent analysis of the low-mass X-ray binary pulsar EXO 0748-676 observational 
data Ozel (2006) concluded that the determination of the mass and radius of this pulsar 
appears to rule out all the soft equations of sate of neutron-star matter. Indeed, multiple 
phenomena have been observed from this low-mass X-ray binary that can be used to deter- 
mine uniquely the mass and radius. 

The three quantities used in Ozel (2006) are the Eddington limit Fem-, the gravitational 
redshift z, and the ratio F cooi /crT c 4 . Following Ozel (2006), we have: 
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Here and in the following we shall adopt natural units where h = c = ks = 1- 

In Eqs. dBJ), ©, and ([5]) G is the gravitational constant, M and R are the stellar mass 
and radius, D the distance to the source, and K e is the electron scattering opacity: 



0.2 (1 + X) cm 2 g- 



(9) 



X being the hydrogen mass fraction of the accreted material. Finally, /oo is the colour 
correction factor which relates the colour temperature T c to the effective temperature T e ff 
of the star. 

Ozel pointed out that the following parameterization of the colour correction factor: 
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leads to an accurate description of the results from model atmosphere calculations (IMadej et al. 2004 ). 

Remarkably, Eqs . (El), (HI), a nd ( [Sj) can be solved to uniquely determine the stellar mass, 
radius, and distance (jOzel 20061 ): 
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The value of the redshift has been reported in Cottam et al. (2002): 

z = 0.35 . 



(11) 
(12) 
(13) 

(14) 



Moreover, the Eddington limit lumi nosity Fem has be en obtained in O zel (2006) by averagin g 
the values determined with RXTE JWolff et al. 2005h and EXO SAT jGottwald et al. 1986h : 

F Edd = 2.25 ±0.23 10~ 8 erg cm- 2 s' 1 . (15) 

Finally, the ratio F coo i/aT^ has been inferred from the RXTE data (jWolff et al. 2005h : 

Fcool 



1.14 ±0.10 (km/kpcY 



(16) 



Using the values in Eqs. (fT4"|) . (|T5l) . and (fTBT) we easily obtain from Eqs. (|TT1) . (|T2^) . and (fT3"|) : 

(17) 



M f 4 
- [0.2009 ± 0.0271] ^ 



M, 



O 



f 4 

[1.315 ± 0.177] km ^ 



P 

D = [1.232 ± 0.137] kpc — 



(19) 



The main uncertainty that affects the mass, radius, and distance of the binary pulsar resides 
in the hydrogen mass fraction X of the accreting material and the colour correction factor 

foo- 

Following Ozel (2006), if we assume for the colour correction factor: 



fo 



1.37 



(20) 




Fig. 1. — Full lines are the gravitational mass M versus the stellar radius R for p-stars 
for different values of \fg~H (in Gev). The dashed lines are the mass-radius curves for 
neutron stars by assuming the stiffest equation of state (MSO) and a softer equation of 
state (AP4, see Lattimer & Prakash 2001). Full triangles correspond to the EXO 0748- 
676 mass-radius constraints assuming ~ 1.37 for two different values of hydrogen mass 
fraction X. Full circles are the best-fitting mass-radius values of the model 2 in Table 1 of 
Poutanen & Gierlinski (2003) for the accreting X-ray millisecond pulsar J1808. The thick 
dashed line corresponds to the mass-radius curve for RXJ 1856 obtained by solving Eq. ( 1271) 
with R°° ~ 6.6 Km. The dotted line is the semi-empirical M-R curve for Her X-l discussed 
in Reynolds et al. (1997) with assumed luminosity of Lx = 3.5 x 10 37 ergs^ 1 . The 
arrows indicate that the M-R curves for Her X-l and RXJ 1856 are to be intended as lower 
limits. The region comprises between the two dot-dahed lines and the M-R curves with 
y/gH = 0.6Gev, 0.35Get> is the allowed region for p-star. 



and using the extreme value of hydrogen mass fraction X ~ 0.70, we get : 

M 



2.10 ± 0.28 , (21) 



R = 13.62 ± 1.84 km . (22) 

On the other hand, for binary systems like EXO 0748-676 a helium-rich companion should 
be expected (see, for instance Lewin 1993). In this case the value X ~ 0.30 seems to be 
more appropriate. If this is the case, we get: 



— = 2.72 ± 0.37 
O 



(23) 
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R = 17.82 ± 2.40 km . (24) 
In Fig. [T]we compare Eqs. (121"]) . (1221). and Eqs. fl23|) . fl24|) with the mass-radius relation ob- 



tained within our p-star theory (ICea 2004al ). We can see that for both hydrogen abundances 
our model is able to account for the inferred values, Eqs. (T2"Tj) . fT2"2"|) and Eqs. (T2"3"j) . (T2"4"l) . to- 
gether with the constraint Eq. (fl4T) . On the contrary, the mass-radius curve for neutron stars 
is consistent with the determination Eq. (12T|) and Eq. (1221) only by assuming the stiffest equa- 
tion of state (labelled MSO in Fig. [fl see Lattimer 2001), which corre sponds to high-densit y 



neutron matter without non- linear vector and isovector interactions (jMuller fc Serot 19961 ). 
Note, however, that even the stiffest high-density neutron matter equation of state is ruled 
out if we consider the more realistic hydrogen mass fraction X ~ 0.30. 

Interestingly enough, we find that regardless of the actual value of hydrogen mass frac- 
tion the distance of the binary system is quite tightly constrained: 

(8.87 ± 0.99) kpc < D < (12.38 ± 1.39) kpc . (25) 

We see that Fig. [T] shows that EXO 0748-670 can be accounted for by our p-star model 
irrespective on the assumed value of the hydrogen mass fraction. 

We see that the stiffest equation of state for neutron-star matter is marginally compatible 
with observations from EXO 0748-670. However, we argue in a moment that there is some 
tension with several pulsar data which would require a softer equation of state. 

It has been proposed long time ago that t he compact ac creting object in the famous 



X-ray binary Herculses X-l is a strange star (ILi et al. 19951 ). This proposal was based 
on the comparison of a phenomenological mass-radius relation for Herculses X-l (see for 
instance Shapiro & Teukolsky 1983) with theoretical M-R curves for neutron and strange 
stars. The analysis in Li et al. (1995) has, however, been criticized by Reynolds et al. 
(1997). These authors, using a new mass estimate together with a revised distance, which 
leads to a somewhat higher X-ray luminosity, argued that the hypothesis that Herculses 
X-l is a neutron star is not disproved. As a matter of fact, Reynolds et al. (1997) found 
that there is marginal consistency with observations if one adopts for neutron stars a very 
soft equation of state. At the same time, these authors pointed out that the hypothesis of a 
strange star can be ruled out since the theoretical curves no longer intercept the observational 
relations within the permitted mass range. However, it is evident from Fig. 1 in Reynolds 
et al. (1997) that the mass-radius constraints for Her X-l seem to be incompatible with 
the stiffest equation of state displayed in our Fig. [1] where we display the semi-empirical 
M-R curve for Her X-l discussed in Reynolds et al. (1997) corresponding to the luminosity 
Lx = 3.5 x 10 37 erg s -1 . On the other hand, the phenomenological mass-radius relation for 
Herculses X-l could be compatible with a softer equation of state like the one labelled AP4 
in Fig. [1] (for the nomenclature on the equations of state see Lattimer & Prakash 2001). 
The millisecond pulsar SAX J1808-3658 (henceforth J1808) is one of the most studied 



- 8- 



accreting pulsar (for instance, see Van der Klis 1995). Recently, Poutanen & Gierlinski 
(2003) studied the pulse profile of J1808 at different energies. In particular, they derived 
simple analytical formulae for the light curves. By fitting the observed pulse profiles in 
the energy band 3 — 4 keV and 12 — 18 keV Poutanen & Gierlinski (2003) constrain the 
compact star mass and radius. The best-fitting parameters for two different models are 
presented in Table 1 of Poutanen & Gierlinsk i (2003). It turns out that t he results from 
the two adopted models are quite consistent (IPoutanen fc Gierlinski 20031 ). In Fig. [I] we 
report the results form model 2 in Poutanen & Gierlinski (2003). Again, we see that the 
mass-radius constraints for J1808 are more consistent with a soft equation of state for high 
density neutron matter. On the other hand, our p-star model can easily account for the 
observed mass-radius values for the accreting X-ray millisecond pulsar J1808. 

Remarkably, our previous conclusion is confirm ed by the recent a nalysis of the light 
curves of J1808 during its 1998 and 2002 outbursts (ILeahy et al. 20071 ). Indeed, Leahy et 
al (2007) obtain that at the 3 a level the radius must satisfy R < 11.9 km and the mass 
M < 1.56 M Q . 

RXJ 1856.5-3754 (in the following RXJ 1856) is the nearest and brightest of a class of 
isolated radio-quiet compact stars. RXJ 1856 has been observed with Chandra and XMM- 
Newton (jBurwitz et al. 20031 ) , showing that the X-ray spectrum is accurately fitted by a 
blackbody law. Assuming that the X-ray thermal emission is due to the surface of the star, 
Burwitz et al. (2003) found for the effective radius and surface temperature: 



R c 



AA 



Km , 

120 pc 



63 eV 



where 



R c 



R 



T 



2GM 
R 



2GM 
R 



(26) 



(27) 



It should be stressed, however, that in the observed spectrum there is also an optical emission 
in excess over the extrapolated X-ray blackbody. By interpreting the optical emission as a 
Rayleigh- Jeans tail of a thermal blackbody emission, one finds that the optical and X data 
can be also fitted by the two blackbody model. In this case, the spectral fit yield an effective 
radius (jBurwitz et al. 20031 ): 



R 



OO ^> 
no 



16 Km 



d 



(28) 



120 pc 

The two blackbody model, however, does not furnish an acceptable description of the ob- 
served spectrum. Indeed, interestingly enough, quite recently the distance measurement o f 
RXJ 1856 has been reassessed and it is now est imated to be at about 180 pc fcaplan 2003h . 
Moreover, from recent parallax measurements (IKaplan 20031 : Ivan Kerkwijk Kaplan 20071 ) 
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we infer that there is a lower limit to the distance of RXJ 1856: 

d > 160 pc . (29) 

In fact, Eq. (J29j) excludes the two-blackbody interpretation, for in this case from Eq. ff28l 
we should obtain: 

R°° > 21 Km , (30) 

which is too large for a neutron star. In addition, if we consider the second nearest 
isolated com pact star RXJ 072 .4-3125 (henceforth R XJ 0720), which has been detected 



by ROSAT (Habcr l et al. 1997: Motch fc Haberl 19981 ) and observed with XMM-Newton 



( [Cropper et al. 20011 : iPaerels et al. 20011 ). then we see that the spectrum of RXJ 0720 is al- 
most identical to that of RXJ 1856. Indeed, RXJ 0720 exhibits a blackbody X-ray spectrum, 
a large X-ray to optical flux ratio, and an optical emission in excess over the extrapolated 
X-ray black body. In this case, however, the recent analysis of the optical, ultraviolet, and 



X-ray data (IKaplan et al. 20031 ) showed that the optical spectrum of RXJ 0720 is not well 
fitted by a Rayleigh- Jeans tail, but it is best fitted by a non thermal power law. This strongly 
suggest that the pulsar surface emission cannot account for the observed soft emission spec- 
tra, and this last emission must be magnetospheric in origin. For instance, i n Cea (2004b) 



we suggested that the soft spectrum originates from synchrotron radiation ([Wallace 19771 ) 
emitted by electrons with power-law energy spectrum. In this case, the radiation in the soft 
spectrum should display a rather large linear polarization. 

Our previous discussion lead to conclude that most realistic interpretation of the X- 
ray spectrum is thermal emission due to the whole surface of the star. Now, assuming 
d ~ 180 pc, from Eq. (1261) we get R°° ~ 6.6 Km. Fro m this value of R°° we can solve 



Eq. (!27j) to constrain the mass and radius of RXJ 1856 (ICea 2004bl ). The result of this 
analysis, displayed in Fig. [TJ indicates that there are stable p-star configurations which agree 
with observational data. However, it should be keep in mind that R°° is a lower limit of 
the true stellar radius. Nevertheless, from Fig. [TJ we see that the mass-radius curve for RXJ 
1856 is more consistent with a soft equation of state. 

Finally, from Fig. [TJ assuming that the observed compact objects are p-stars we infer 
that the chromagnetic condensate is constrained in the rather narrow interval: 



0.35 Get; < xj gH < 0.6 Gev . (31) 

As a consequence, the allowed region for p-stars is the region in Fig. [TJ bounded by the two 
dot-dahed lines and the M-R curves with y/gH = 0.6Gev, 0.35Gev . 
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3. COOLING 

The thermal evolution of compact stars is an important tool to investigate the state 
of dense matter at supra- nuclear densities. In fact, observations of the thermal photon 
flux emitted from the surface of the stars provide valuable information about the physical 
processes operating in the interior of these objects. 

The cooling in p-stars has been discussed for the first time in Cea (2004a). As for neutron 
stars, the predominant cooling mechanism of newly formed p-stars is neutrino emission. In 
p-stars neutrino cooling dominates for about 10 2 — 10 3 years. Subsequently, after the internal 
temperature has sufficiently dropped, the photon emission ov ertakes neutrinos. 



Let us briefly review the cooling in p-stars (ICea 2004al ). Assuming stars of uniform 
density and isothermal, the cooling equation is: 

dT 

C v — = - (L u + L 7 ) , (32) 

where L v is the neutrino luminosity, L 7 is the photon luminosity and Cy is the specific heat. 
Assuming blackbody photon emission from the surface at an effective surface temperature 
T s we get: 

L 7 = 4ttR 2 ctsbTI, (33) 

where <jsb is the Stefan — Boltzmann constant. Following Shapiro & Teukolsky (1983), in 
Cea (2004a) we assumed that the surface and interior temperature were related by: 

^ = 10- 2 a , 0.1 < a < 1.0 . (34) 
Eq. (134")) is relevant for a p-s tar which i s not bare, namely for p-stars which are endowed 



with a thin crust. It results (ICea 2004bl ) that p-stars have a sharp edge of thickness of the 
order of about 1 fermi. On the other hand, electrons which are bound by the coulomb 
attraction, extend several hundred fermis beyond the edge. It follows, then, that on the 
surface of the star there is a positively charged layer which is able to support a thin crust 
of ordinary matter. The vacuum gap between the core and the crust of order of several 
hundred fermis leads to a strong suppression of the surface temperature with respect to the 
core temperature. In principle, the actual relation between Ts and T can be obtained by 
studying the crust and core thermal interactions. In any case, our phenomenological relation 
Eq. (1341) allows a wide variation of Ts, which encompasses the neutron star relation (see, for 
instance, Gudmundsson 1983). Moreover, as we discuss below, our cooling curves display a 
rather weak dependence on the parameter a in Eq. (1341) . 

In Cea (2004a) we showed that the dominant cooling processe s by neutrino emission are 



the direct /3-decay quark reactions (jlwamoto 19801 : iBurrows 19801 ): 



d — > u + e + v e , u + e — > d + v e . (35) 
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We find the following neutrino luminosity (jCea 2004al ): 



where Tg is the temperature in units of 10 9 °K, and 6q — 2.51 10 14 g cm~ 3 is the nuclear 
density. So that, for typical values of parameters we have: 

L v ~ 10 36 ^ T 9 8 . (37) 

Note that the neutrino luminosity L v has the same temperature dependence as the neutrino 
luminosity by modified URCA reactions in neutron stars (see, for instance Shapiro & Teukol- 
sky 1983), but it is more than two order of magnitude smaller. This peculiar dependence on 
the core temperature is due to the presence of the strong chromomagnetic condensate which 
strongly constraints the quark transverse motion. 

It is worthwhile to stress that, since our neutrino luminosity is reduced by more than 
two order of magnitude with respect to neutron stars, the maximum allowed quiescent lumi- 
nosity of isolated pulsars is about tw o order of magnit ude greater than the maximum allowed 



surface luminosity in neutron stars (IVan Riper 19911 ). Thus, while our theor y allows to ac 



count for luminosities up to 10 36 erg s _1 as observed in magnetars (jCea 20061 ). the standard 
model based on neutron stars is in striking contradiction with observations. 

A further support on our neutrino luminosity Eq. (1371) comes from superbursts, namely 
rare, extremely energetic, and long duration type I X-ray bursts, from low mass X-ray bina- 
ries (see Kuulkers 2004). Indeed, it has been pointed out that to account for the observed 
i gnition depths in superbursts one needs slow ne utrino emission processes with emissivity 



(IPage & Cumming 20051 : ISteiner & Madsen 2006T ): 



e„ = Q v If -^f- , Q v ~ 10 18 - 10 22 . (38) 
cm 4 s 

It is evident that the phenomenological neutrino emissivity Eq. (138|) needed to explain ob- 
served superbursts from low mass X-ray binaries cannot be accounted for within the standard 
neutron star model. On the contrary, from our Eq. f)37p we infer for the neutrino emissivity 
in p-stars: 

e v ~ 10 18 Tg , (39) 

cm 6 s 

which naturally fits in the phenomenological allowed range of values. 



Finally, to determine the thermal evolution of p-stars we need the specific heat (jCea 2004al ) 



Oy * 0.92 10- T, M_ £! ( ^L? , (40) 
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Fig. 2. — Comparison of p-star cooling curves with several data for the effective surface 
temperature. Different lines correspond to compact star mass values indicated in the legend 
and different values of the parameter a in Eq. (jMj) . Data points with error bars are taken 
from Yakovlev & Pethick (2004) and references therein. 



which, for typical parameter values, in physical units reads: 

C v ~ 10 39 ^ T 9 . (41) 
Note that, from Eq. (j4Tl) i t follows that the p-star sp ecific heat is of the same order of the 



neutron star specific heat (IShapiro fc Teukolsky 19831 ). To obtain the thermal history of a 
p-star we integrate Eq. (|32|) by assuming the initial temperature = 1.4 ( Cea 2004a ). 



We stress, however, that the thermal history is almost independent on the assumed initial 
temperature as long as T^' < 10 10 °K. 

In Fig. [2] we report our cooling curves for three different values of the parameter a in 
Eq. fl34l) . It is worthwhile to note that the cooling curves are almost independent on the 
p-star mass. Moreover, there is a weak dependence on a up to age r ~ 10 3 years. We 
compare our theoretical cooling curves with several pulsar data taken from the literature. 

In any case we see that the agreement between theoretical cooling curves and observa- 
tional data is quite satisfying. In particular, we see that, at variance with neutron stars, 
our peculiar neutrino luminosity allows sizeable surface effective temperature, up to 10 5 °K, 
for compact stars with age r > 10 6 years. Indeed, Fig. [2] suggests that this distinguishable 
feature of the p-star model is corroborated by observations. 
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4. CONCLUSIONS 



The proposal for p-stars originates from our non-perturbative investigations of QCD 
(Cea & Cosmai 2003, 2005) which suggested that quarks and gluons get deconfined in strong 
enough chromomagnetic fields. This, in turns, leads us to argue that the deconfined QCD vac- 
uum is characterised by long-range chromomagnetic correlations and that p-matter, namely 
almost massless up and down quarks immersed in a chromomagnetic condensate, is formed 
in the collapse of the core of an evolved massive star (ICea 2004al ). In addition, we already 
argued that p-stars are more stable than neutron stars whatever the value of the chromomag- 
netic condensate. As a consequence, the true ground state of QCD in strong gravitational 
fields is not hadronic matter, but p-matter. In other words, the final collapse of massive 
stars leads inevitably to the formation of p-stars. 

The results discussed in the present paper, indeed, indicate that p-stars are able to re- 
produce in a natural manner several observations from pulsars. On the other hand, we have 
seen that within the standard model based on neutron stars there is some tension in the 
determination of the equation of state for high density neutron matter. However, we would 
stress that mass-radius constraints are affected by important systematic errors and are not 
therefore so constraining. Thus, we need more precise mass-radius determinations to reach 
a firm conclusion. 
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